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Summary. — We study the Wigner functions of the nucleon which provide multi- 
dimensional images of the quark distributions in phase space and combine in a single 
picture all the information contained in the generalized parton distributions (GPDs) 
and the transverse-momentum dependent parton distributions (TMDs). In particu- 
lar, we present results for the distribution of unpolarized quarks in a longitudinally 
polarized nucleon obtained in a light-cone constituent quark model. Treating the 
Wigner distribution as it was a classical distribution, we obtain the quark orbital 
angular momentum and compare it with alternative definitions given in terms of the 
GPDs and the TMDs. 

PACS 12.38.-t - Quantum chromodynamics. 
PACS 12.39.-x - Phenomenological quark models. 
PACS 14 . 20 . Dh - Protons and neutrons. 



1. Introduction 

Quark Wigner distributions provide joint position-and-momentum (or phase-space) 
distributions, encoding in a unified picture the information obtained from transverse- 
momentum dependent parton distributions (TMDs) and generalized parton distribu- 
tions (GPDs) in impact-parameter space. The concept of Wigner distributions in QCD 
for quarks and gluons was first explored in refs. [1, 2}. Neglecting rclativistic effects, the 
authors introduced six-dimensional (three position and three momentum coordinates) 
Wigner distributions. In a recent work [3], we used the connection between Wigner dis- 
tributions and generalized transverse-momentum dependent parton distributions (GT- 
MDs) [4] to study five-dimensional distributions (two position and three momentum 
coordinates) in the infinite-momentum frame providing a picture consistent with special 
relativity. However, it is well known that the quantum-mechanical phase-space distribu- 
tions do not have a density interpretation because of the uncertainty principle. Accord- 
ingly, Wigner distributions are not positively defined. Nevertheless, the physics of the 
Wigner distributions is very rich and one can select certain situations where a semiclas- 
sical interpretation is still possible. 



© Societa Italiana di Fisica 



1 



2 



C. LORCE and B. PASQUINI 



The purpose of this contribution is to investigate the phenomenology of the quark 
Wigncr distributions. As a matter of fact, since it is not known how to access these 
distributions directly from experiments, phenomcnological models are very powerful in 
this context. Collecting the information that one can learn from quark models which were 
built up on the basis of available experimental information on GPDs and TMDs, one can 
hope to reconstruct a faithful description of the physics of the Wigner distributions. To 
this aim we will rely on models for the light-cone wave functions which have already been 
used for the description of the GPDs [5, 6], the TMDs [6, 7, 8, 9, 10] and electroweak 
properties of the nucleon [6, 11, 12, 13, 14]. 

2. Wigner distributions 

Similarly to refs. [1, 2], we define the Hcrmitian Wigner operators for quarks at a 
fixed light-cone time y + = as follows 

(1) W^{% u %^x) = \ J dZ (2 ^ e^ +z ~- fc "^> ff(g - f )TW^(y + f )\ z+=Q , 

with j/ M = [0, 0, 6j_], p + the average nucleon longitudinal momentum and x = k + /p + 
the average fraction of nucleon longitudinal momentum carried by the active quark. 
The superscript T stands for any twist-two Dirac operator T = 7 + , 7 + 75, i<J J+r )§ with 
j = 1, 2. A Wilson line W = W(y— |,j/ + ||n) ensures the color gauge invariance of the 
Wigncr operator. In the following we will focus on the quark contribution, ignoring the 
gauge-field degrees of freedom and therefore reducing the gauge link W to the identity. 

We define the Wigner distributions p' r '(6j_, k±,x, S) in terms of the matrix elements 
of the Wigner operators (1) sandwiched between nucleon states with polarization S 

(2) pW(b ± ,k ± ,x,8) = J^±( p +,^,S\W^(b ± ,k ± ,x)\p+,-^J). 

As outlined in rcf. [3] , such matrix elements can be interpreted as two-dimensional Fourier 
transforms of the GTMDs in the impact-parameter space. Although the GTMDs are in 
general complex-valued functions, their two-dimensional Fourier transforms are always 
real-valued functions, in accordance with their interpretation as phase-space distribu- 
tions. We note that, like in the usual quantum-mechanical Wigner distributions, b± 
and k± are not Fourier conjugate variables. However, they are subjected to Hcisen- 
berg's uncertainty principle because the corresponding quantum-mechanical operators 

do not commute [b±,k±] ^ 0. As a consequence, the Wigner functions can not have a 
strict probabilistic interpretation. There are in total 16 Wigner functions at twist-two 
level, corresponding to all the 16 possible configurations of nucleon and quark polar- 
izations. Here we will discuss only one particular case, namely the distortion in the 
distribution of unpolarizcd quarks due to the longitudinal polarization of the nucleon 
Plu — /°' 7 {b±, k±, x, +e z ) — p' 7 '(b±, k±, x, —e z ) which has a close connection with the 
quark orbital angular momentum (OAM). Other configurations for the quark and nucleon 
polarizations can be found in ref. [3]. 

In fig. 1, the upper panels show the distortions in impact-parameter space for u 
(left panels) and d (right panels) quarks with fixed transverse momentum kj_ = k± e y 
and k± = 0.3 GcV, while the lower panels give the corresponding distortions in the 
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Fig. 1. - The distortions of the x- integrated Wigner distributions of unpolarized quarks due to 
the spin of the proton (pointing out of the plane) in a LCCQM. Upper panels: distortions in 
impact-parameter space with fixed transverse momentum k± = k± e y and k± = 0.3 GeV. Lower 
panels: distortions in transverse-momentum space with fixed impact parameter b± — b± e y and 
b± = 0.4 fm. The left (right) panels show the results for u (d) quarks. 



transverse-momentum space with fixed impact parameter b± = b± e y and b± = 0.4 fm. 
We observe a clear dipole structure in both these distributions which indicates that the 
(quasi-)probability for finding the quark orbiting clockwise is not the same as for the 
quark orbiting anticlockwise, leading in average to a nonvanishing OAM. We learn from 
these figures that the OAM of u quarks tends be aligned with the nucleon spin, while the 
OAM of d quarks tends be antialigned with the nucleon spin. In particular, we notice 
that the distortion induced by the quark OAM is stronger in the central region of the 
phase space (k± <C and b± ^;), for both u and d quarks. The distortion in the 6j_-space 
(see upper panels of fig. 1) is more extended for d quarks than for u quarks, whereas the 
opposite behavior is found for the distortion in the fc^-space (see lower panels of fig. 1). 
In the case of d quarks, we also observe a sign change of the distributions in the outer 
regions of phase space (k± ^> and b± 3>) which corresponds to a flip of the local net 
quark OAM. Note that such a direct access to quark OAM is not possible with GPDs and 
TMDs since none of them describe at leading twist the distortion in the distribution of 
unpolarized quarks due to the longitudinal polarization of the nucleon. This is because 
one needs the correlation between b± and k± which is lost by integrating over b± or k± . 
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3. Quark orbital angular momentum 

If one neglects gauge-field degrees of freedom, the quark OAM operator for a given 
flavor q can be unambiguously defined as 

* 



(3) L9z = ~2j d r ^ 7 [fxVr)^\ 

where normal ordering is understood. We define the quark OAM t\ as the following 
matrix clement of the quark OAM operator 

where the incoming and outgoing momenta of the nucleon are given by p = P — ^A 
and p' = P + iA, respectively. Since the nucleon state |p, A) in eq. (4) is normalized as 
(p', +|p, +) = 2P+ S(A+) (2tt) 3 <5( 2 )(Aj_), the forward limit p' = p has to be treated with 
care. Integrating over A + and Aj_ leads to consistent expressions while avoiding the 
use of normalizable wave packets or infinite normalization factors. It turns out that this 
definition of quark OAM coincides with treating the Wigner functions as if they were 
classical distributions, and obtaining the quark OAM by calculating the integral over 
phase space of the quark distribution in a longitudinally polarized nucleon multiplied by 
the naive OAM operator (b±_ x k±) z [15] 



(5) l\ = fdxd 3 k±d 2 b± (b± x fcj_) p^ q (b ± ,k A 



3'1. Overlap representation. - Following the lines of refs. [16, 17], we obtain an overlap 
representation of eq. (5) in terms of light-cone wave functions (LCWFs). Since the 
quark OAM operator is diagonal in light-cone helicity, flavor and color spaces, we write 
the quark OAM as the sum of the contributions from the A-parton Fock states t\ = 



N N 



(6) l?"* = \\ [dx] N [d 2 kA_] N f2^ E ~ *«) [*^» Vfe„) Kp(r) 



Recently, it has been suggested, on the basis of some quark-model calculations, that 
the transverse-momentum dependent distribution h\ T may also be related to the quark 
OAM [18, 19]. It is however important to note that in this case, the quark OAM is 
defined in more naive terms 

. - N 

(7) C™« = -\ j [dx] N [d^] N ^ 

Finally, according to Ji's sum rule [20], the quark OAM can also be extracted from 
the following combination of GPDs 

(8) L ^ = ^/ dx { x l H9 ( x ' > ) + E9 ( x >°> )]- ^ 9 (a;,0,0)} . 
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The overlap representation of the ./V-parton Fock state contribution is then given by 

i r N 

(9) LW>* = - / [dx] N [d 2 k ± ] N J2Sm 
- 1 i=i 

\ (xi - Xi) *^(r) 2 + Mxi (A 



3'2. Model results and discussion. - We present in Table I the results from a LCCQM 
and the light-cone version of the chiral quark-soliton model (LCxQSM) restricted to the 
three-quark sector [3]. 

Table I. - The results for quark orbital angular momentum (see eqs. (6), (7) and (9)) from a 
LCCQM and the LCxQSM for u-, d- and total (u + d) quark contributions. 



Model 

q 


u 


LCCQM 
d 


Total 


u 


XQSM 
d 


Total 


£ q 


eq. (6) 


0.131 


-0.005 


0.126 


0.073 


-0.004 


0.069 




eq. (7) 


0.169 


-0.042 


0.126 


0.093 


-0.023 


0.069 


LI 


eq. (9) 


0.071 


0.055 


0.126 


-0.008 


0.077 


0.069 



First of all, we see that all definitions give the same values for the total quark OAM. 
This is expected since once summed over all flavors the expressions (6), (7) and 
(9) become identical. We see also that there is more net quark OAM in the LCCQM 
(Eg 1 '! = °- 126 ) than m the LCxQSM (Eg^I = 0.069). However, all definitions dis- 
agree for the separate flavor contributions. Both the LCCQM and the LCxQSM predict 
that l\ (from the Wigner functions) and C% (from the h^ T TMD) are positive for u 
quarks and negative for d quarks, with the u-quark contribution larger than the d-quark 
contribution in absolute value. For L q z (from the Ji's sum rule) the LCCQM predicts 
the same positive sign for the u and d contributions, with the isovector combination 
L" — L d z > 0, similarly to a variety of relativistic quark model calculations. Instead, the 
LCxQSM gives L z < and L z > 0, and thereefore L z — L z < 0, in agreement with 
lattice calculations [21]. 

It is surprising that l\ =^ L q z since it is generally believed that Jaffe-Manohar and 
Ji's OAM should coincide in absence of gauge degrees of freedom [22]. Note that a 
similar observation has also been made in the instant- form version of the xQSM [23]. 
One may argue that this comes from a problem on the model side (e.g. spurious surface 
contributions or violation of rotational symmetry). This is unlikely considering that the 
isoscalar part is the same l z = L z . It is only for the separate contributions that the 
different definitions disagree. Unless rotational symmetry implies non-trivial relations 
between the different components of the LCWFs, one should expect in general different 
results for (6) and (9). We emphasize also that the different contributions to OAM 
depend on the reference point with respect to which OAM is defined. While it is clear 
that £ q z and C q z are defined with respect to the transverse center of momentum and origin 
of axes, respectively, it is not clear from the overlap representation (9) with respect to 
which point Ji's OAM is defined. 
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4. Conclusion 

In summary, we have presented the first model calculation of the Wigner distributions 
within the light-cone formalism, using a derivation which is not spoiled by relativistic 
corrections. The results within a light-cone constituent quark model and the light-cone 
version of the chiral quark-soliton model are very similar, and allowed us to sketch some 
general features about the behavior of the quarks in the nucleon when observed in the 
6_i_-plane at fixed k±, or in the fcj_-plane at fixed b±. We discussed also three different 
definitions of quark OAM and compared them using the above two models. While all 
definitions agree on the isoscalar combination, they disagree on the isovector one. Possible 
reasons for this observation have been discussed. 

* * * 
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